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We propose a new solution to the origin of dark energy. We suggest that it was created dynamically from the 
condensate of a singlet neutrino at a late epoch of the early Universe through its effective self interaction. This 
singlet neutrino is also the Dirac partner of one of the three observed neutrinos, hence dark energy is related to 
neutrino mass. The onset of this condensate formation in the early Universe is also related to matter density and 
offers an explanation of the coincidence problem of why dark energy (70%) and total matter (30%) are comparable 
at the present time. We demonstrate this idea in a model of neutrino mass with (right-handed) singlet neutrinos 
and a singlet scalar. 



The astrophysical observations that the bary- 
onic and dark matter together account for only 
about 30% of the total matter while the remain- 
ing 70% have the property of producing negative 
pressure, interpreted as the dark energy, remains 
a challenging problem in cosmology. A possible 
explanation that the dark energy is simply the 
vacuum energy given by the standard model scale 
fails by many orders of magnitude. It also can 
not explain the fact that the dark energy has been 
comparable to the density of ordinary matter dur- 
ing the recent epoch in the evolution of the uni- 
verse. While solutions to this problem have been 
provided, notably through the presence of a scalar 
held called quintessence[l],[2]the observation that 
the magnitude of the dark energy is comparable 
to neutrino masses suggests to us that the ex- 
planation of the dark energy puzzle may possibly 
lie in the direction of neutrinos, while the other 
large effects coming from electroweak breaking 
and QCD are canceled out by some unknown dy- 
namical effects It is this avenue which we follow 
in our discussion below. 

Neutrinos have been invoked in the past in the 
dark energy problems. For example, recently it 
has been argued in a theory with scalar fields, 



called accelerons, that the dark energy can be 
obtained dynamically if the neutrino masses are 
allowed to vary with the scalar field [3], [4]. We 
propose, however, a somewhat more direct con- 
nection through the formation of neutrino con- 
densates [5], [6], [7], [8]. We elaborate below 

Our starting point is an extension of the stan- 
dard model which includes right-handed neutri- 
nos ViR, i — 1,2,3, in which the neutrino masses 
have the seesaw structure [9] with both the Dirac 
and Majorana masses of the order of eV. In this 
model we introduce a light real singlet scalar <!?, 
so that the effective Lagrangian after the elec- 
troweak symmetry breaking is given by 

C = m Dta iiV aR + M a V aR V aR + fafjVaRVpRS (1) 

where m_Di Q is the Dirac mass. The Majorana 
masses of the right-handed neutrinos M a are as- 
sumed to be real and diagonal which is achieved 
without loss of generality through the choice of 
our basis. This Lagrangian has a discrete sym- 
metry, (— 1) L , where L is the lepton number, so 
that ti and v aR are odd under this symmetry. 

We take the Dirac neutrino masses to be of the 
order of 0.1 eV, while the Majorana masses of the 
right-handed, singlet, neutrinos are assumed to 
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be of the order of 0.001 eV giving rise to pseudo- 
Dirac neutrinos with very small mass differences. 
In order to be consistent with the solar neutrino 
data, which does not allow very small mass dif- 
ferences[10], however, two of the neutrinos are al- 
lowed to be only either Dirac or Majorana, but 
not pseudo-Dirac. The Dirac case will correspond 
to three left-handed and three right-handed neu- 
trinos with the mass matrix given by 
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The two physical neutrino eigenstates {u\ and v-i) 
are the Dirac neutrinos responsible for the solar 
neutrino oscillations, and only the state ^3 is a 
pseudo-Dirac neutrino. 

We consider here the second, Majorana, pos- 
sibility where there is only one sterile neutrino, 
which we assume forms a condensate. The corre- 
sponding neutrino mass matrix will be 



M v = 



/mi \ 

m 2 * 

m 3 

\ m 3 M J 



Here v\ and v 2 are Majorana neutrinos, needed to 
explain the solar neutrino data, while v% is now a 
pseudo-Dirac neutrino. For the rest of our anal- 
ysis, we will restrict our discussions to only the 
third eigenstate, the pseudo-Dirac neutrino with 
a mass matrix 



M v 



777,3 
7773 M 



and represent the left-handed and right-handed 
states as vi and vr. Below we outline the dy- 
namics behind the condensate formation which 
we argue is the source for the dark energy. 

The exchange of the real scalar, S, is found to 
give rise to attraction between the right handed 



neutrinos vr. Thus, in the context of cosmic evo- 
lution, as the universe cools down to a temper- 
ature below the mass of the neutrinos, this at- 
tractive interaction then causes the right-handed 
neutrinos to form condensates, the candidates for 
the dark energy. 

The four-Fermion effective contact interaction 
generated through S exchange is given by 



Hi = -C {y M v M ) (v m v m ) ■ 



(3) 



where the right-handed Majorana neutrinos vm 
have been defined as 



VM = ^R + V C L , 



(4) 



where v c a L is the CP conjugate of v a R and A 
is the Majorana phase, so that the Majorana 
neutrino satisfies the Majorana condition v c M = 
\*vm .A number operator for Majorana particles 
can not be defined, but in cosmology one can de- 
fine the number density per comoving volume of 
any Majorana particle, when the interaction rate 
of the particle is slower than the expansion rate of 
the universe. This allows us to define the chemical 
potential for the particles. Thus, for the neutrino 
condensation formation the decoupling tempera- 
ture may be considered as the cut off scale for the 
neutrinos. Below this temperature, the four-fermi 
form will be valid with the coupling strength C 
given by 



C = 



(5) 



where 7775 is the mass of the scalar field and the 
generation index of the coupling constant / has 
been suppressed . 

Below we describe the formalism for the vr 
condensation formation, ft is convenient to work 
in the Weyl representation, in which the 75 matrix 
is diagonal. The left-handed and right-handed 
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neutrinos can be written as 

o 



v L = (1 - 7 5 ) v 



vr = (1 + 7 ) v 
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"'i = "// (6) 
so that 

vm = 

The four-fermion Hamiltonian can then be ex- 
pressed in terms of the component fields as 

Hi = -C A* 2 xl Xa x\ Xb + xl Xa x\ Xb 

+ xl Xa xl Xb + A 2 xl Xa Xb Xb ] > ( 7 ) 

and the condensate will correspond to a spin-0 
pairing: 

(Xa xl) =e ab D, (8) 

giving us the interaction Hamiltonian in the 
mean-field approximation as: 



H MF = _ 2C 



A* 2 xl Xb D + A 2 xl X h D* 



tab ■ 

(9) 



In terms of the creation and annihilation opera- 
tors of any Majorana field: 

Mx) = J2 \/lf (Usu ps e- ipx + \*flv ps e^) , 

p.s 

(10) 

where e 2 = p 2 + m 2 a is square of the energy of the 
physical pseudo-Dirac neutrinos. The interaction 
Hamiltonian can now be written as: 
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(11) 



The complete Hamiltonian will be a sum of the 
free Hamiltonian (Hq) and the interaction Hamil- 
tonian {Hf IF ).lt can be transformed to a canon- 
ical form 



n = Y, E { b l^+ b li b Pi) > 



(12) 



by a time-dependent transformation. Following 
the standard condensed matter formalism one 
finds that a consistent solution, that allows a non- 
vanishing condensate D^O, gives 
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where k is related to the gap parameter, A and 
[i is the chemical potential. The lower limit on 
the energy integral is taken to be (M a — fi), and 
the upper limit is the cut off scale A, which is 
the decoupling scale for the neutrinos. This limit 
is determined by the requirement that above this 
temperature neutrino interactions are in equilib- 
rium and the number density and the chemical 
potential are not defined. 

The solution of (13) gives us the magnitude of 
the gap: 



A = 2 



2A 



1/3 



(14) 



where x = 27r 2 /[CM Q (37r 2 n I ,) 1/3 ] and the critical 
temperature and the Pippard coherent length are 
given by 
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(15) 
(16) 



In the present example, the coupling of the neu- 
trinos to S becomes strong in the non-relativistic 
limit, which will imply that the condition M a <C 
ms is required to be satisfied for equation (3) 
to remain valid. However, this condition may 
be relaxed considerably in a relativistic treat- 
ment of superconductivity. A numerical study 
of fermions interacting with a scalar field in a 
strong coupling regime i.e. / ~ 1, shows that 
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the scalar-field propagator can be scaled as 1 / m| 
even for M a > mg [11]. In this strong cou- 
pling regime, the right-handed neutrino conden- 
sates thus formed start dominating the universe 
when the size of Cooper pair becomes compara- 
ble to the inter-particle spacing. For £ = 0.1 cm 
- (2 x lO^eV)- 1 and n v ~ 110, we get x ~ 13.5, 
m s ~ 4.6 x 10" 4 eV and A - 4 x 10~ 5 eV. The 
existence of a finite, non zero, gap provides the 
evidence for a condensate 

This condensate, we would like to argue, is a 
dark energy candidate. Since the Cooper pairs 
are formed around the scale of neutrino masses, 
the amount of dark energy density becomes com- 
parable to the matter density in this scenario. 
The amount of dark energy is determined by 
the Major ana mass of the neutrinos, which is 
M ~ 10 -3 eV, and we get the correct order of 
magnitude. Thus, without invoking any dynami- 
cal field like the quintessence or accelerons, we 
have found a natural explanation as to why the 
scales of dark energy and neutrino masses are 
comparable and why dark energy dominates in 
this epoch. Furthermore, our model is also consis- 
tent with the fact that the amount of dark energy 
is the same as the matter density. 

We conclude with a short discussion of the ba- 
sic dynamics of the condensates, which we call 
£.5. The Lagrangian for £ s is given by 

c = (do + in s ) it (do - i^) ^ - v&), 

(17) 

where m and p s represent the mass and the chem- 
ical potential of the condensate respectively, and 

ne s )-m 2 |6| 2 + 5 |6| 4 . 

For the present case we assume, m ~ 2m„ and 
p s — P- We note that the interaction lagrangian 
can be written as 



V(C s )=ml\Ut\+C(v c R VR + hc), 



(18) 



1 It is interesting to note what happens if we consider the 
neutrinos in isolation and discuss their bound state in a 
non-relativistic framework. Because of the miniscule mass 
of the scalar, S, one can approximate the interaction to 
be Coulomb-like . Since / ~ 1, the Bohr radius of the 
resulting atom will be ~ and the binding energy 

~ M a which is essentially of the order of the dark energy 



which gives V(£ s ) after taking into account the 
quartic self-interactions. 

In absence of chemical potential, the equation 
of state for the scalar-field becomes 



P 
P 



KE + V(£ S ) 



For KE <C V(i s ), the scalar field then behaves 
as dark energy with the desired value u> *~ — 1. 
This would imply that g|£ s | 2 <0org<0. The 
coupling constant g can be calculated from the 
knowledge of scattering between the condensates 
from an analogy with atomic physics. One can 
write q = for an attractive interaction with 
a < 0. Under this condition the condensates can 
represent the dark energy [12]. Another possibil- 
ity is the case when there is a finite chemical po- 
tential and there is spontaneous symmetry break- 
ing [8, 13]. In this case one can use the arguments 
of Ref. [8] to estimate the parameters m and g as 
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8^2" 



and g 



N f C 



(19) 



8tt 2 VM" 

In this case the equation of state can be written as 
p = 3p + 4(1 - ^y/l + >/p)y/p [14] which gives 
lu ~ — 1 if p = (^2 — I) 2 and p? > m 2 which 
is consistent with the condensate formation con- 
dition [13]. p and p are made dimensionless by 
factoring out with to 4 /Ag. 

In summary we have proposed a new solution of 
the dark energy problem where the dark energy is 
the condensate formed by self interaction of right 
handed (singlet) neutrinos generated through the 
exchange of a singlet scalar. Since this neutrino 
is the Dirac partner of one of the three observed 
neutrinos, the dark energy is related to the neu- 
trino mass. The fact that the matter and dark 
energy are comparable follows naturally from our 
model. 
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